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ABSTRACT 

An anisotropic creep model is formulated for metallic composites with 
strong fibers and low to moderate fiber volume percent (< 40%). The ideal- 
ization admits no creep in the local fiber direction and assumes equal creep 
strength in longitudinal and transverse shear. Identification of the matrix 
behavior with that of the isotropic limit of the theory permits characterization 
of the composite through uniaxial creep tests on the matrix material. 
Constant and step-wise creep tests are required as a data base. The model 
provides an upper bound on the transverse creep strength of a composite 
having strong fibers embedded in a particular matrix material. Comparison of 
the measured transverse strength with the upper bound gives an assessment of 
the integrity of the composite. Application is made to a Kanthal composite, a 
model high-temperature composite system. Predictions are made of the creep 
response of fiber reinforced Kanthal tubes under interior pressure. 


INTRODUCTION 

Monolithic structural alloys used in high temperature environments 
(T/Tm > .3 or .4) exhibit complex thermomechanical behavior that is both 
time and history dependent. The wealth of research activity over the past 
decade directed toward high temperature structural alloys is testimony to these 
behavioral complexities and to the need for mathematical representations 
(constitutive equations) of this behavior in support of component design. 

When two metallic alloys are combined to form a composite that is 
exposed to high temperature, one or both of the constituents may behave 
inelastically in this same time dependent, hereditary way. Although inelastic 
response may be suppressed under some special conditions, viz., when a single 
stress component is aligned directly with strong fibers, high temperature 
response under general states of stress inevitably involves inelasticity (creep, 



relaxation, rate sensitive plasticity, etc.) Micromechanics, successful in 
predicting composite behavior with elastic constituents, has a difficult task of 
predicting the overall deformation behavior of a composite with its constituents 
behaving in a time dependent, hereditary way, particularly when an interphase 
material develops and its properties are also time and history dependent. 

An alternative is to consider the composite as a material in its own 

right, with its own properties that can be measured for the composite as a 
whole. This amounts to idealizing the composite as a psuedohomogeneous, 
anisotropic material and applying the principles of continuum mechanics. This 
is the approach taken in [1,2,3] and more recently in [4,5], dealing with the 
viscoplastic behavior of unidirectional (one fiber family) metallic composites. 

In this paper, we focus on the creep response of metallic composites and 

specialize some of the results in [5] for the case of very strong fibers and low 

to moderate fiber volume percent (i.e., less than 40% f.v.p). The former 

condition limits the deformation along the fiber direction to be essentially 
elastic. Experiments on a W/Kanthal (~ 35% f.v.p) composite reveal creep 
rates along and transverse to the fiber direction that are different by several 
orders of magnitude [6]. The latter condition ensures domination of inelastic 
response by the matrix and supports an idealization of equal creep strength in 
transverse and longitudinal shear. 

A key feature of the simple theory presented here is that all the 
properties necessary for characterizing creep of a metallic composite can be 
determined through uniaxial creep experiments on the matrix material. Once 
these properties have been established, predictions of creep under variable and 
multi axial stress can be made for the composite. The theory is independent of 
fiber properties requiring only that the fibers are much stronger than the 
matrix material. 

We first formulate the specialized creep theory; then, we outline a 
procedure for determining the material constants from tests on the matrix 
material. This procedure is subsequently applied to a Kanthal composite 
material. Using the constants determined for the Kanthal matrix, predictions 

of creep response of the Kanthal composite are made under transverse tension 
and under biaxial stress states corresponding to reinforced thin walled tubes 
subjected to an interior pressure. 
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FORMULATION OF THE CREEP THEORY 

Earlier work [5] introduced the invariant combination 


♦ = r i + h h 


+ 


V 


4(4o/ 2 -1) ^ 


( 1 ) 


in representing the viscoplastic response of a metallic composite idealized as 
locally transversely isotropic. The invariants Ip ^ and Ig have physical 
meanings, cf. [5], and are defined as 


in which 


h ~ J 2 “ 1 + \ *3 

X 2 - * - h 


( 2 ) 
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and Ey = Sij - ay 


( 3 ) 


where: 

s.j = the components of deviatoric applied stress, a-j = those of the deviatoric 
internal (back) stress and S- • = the components 

J 


of effective stress. d- 


1 


represents a unit vector denoting the local fiber direction (direction of 
transverse isotropy). tj in eq. (1) is the ratio of (creep) strength under 
longitudinal shear stress to that under transverse shear Krp. i.e., 


77 = 



( 4 ) 


Similarly, u is the ratio of strength in longitudinal tension 
transverse tension Yrp, i.e. 



to that in 

( 5 ) 
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As we wish to consider the case of very strong fibers, we take 


OJ - * CD (6) 

Also, reflecting low to moderate fiber volume percent, we take = Krp = K 
so that 

V = 1 (7) 

With eqs. (6) and (7) the invariant (|> becomes 

♦ = J 2 - \ I 2 (8) 

Guided by [5] we propose a creep representation based on the invariant <J> 
and an analogous one defined as 

* = 4 - ! s 2 (») 

where 

4 = S a ij a ji ’ J= D ij a ji ( 10 ) 

A full statement of the proposed (isothermal) theory is as follows 


4 4 
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r ij = V 2 [ ( 3D iJ - *ij ) 
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The functions and constants to be specified in characterizing a particular 
material are: 

f(F), h(G), r(G), a Q and e Q (16) 

Under longitudinal stress (fig. la) this theory predicts no inelastic strain 

rate (no creep), i.e., = a^ = 0. 

Under transverse stress (fig. lb) the theory gives 


F = 


_ _ 

2 


a-a 

, G = 

a 

a o 

a o 


a 


= MG) 


f = m 


-r(G) = h(G)f(F) r(G) 


(17) 

(18) 

(19) 


in which a = a ^ = the transverse applied stress, a = a ^ = the transverse 

component of the internal stress and e e ^2 = the transverse inelastic strain 
rate. 

If, as in earlier work [4, 5], we take 

f(F) = F n , h(G) = H/G^ and r(G) = RG m ~^ 
eqs. (18 k 19) become 


( 20 ) 
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where 


H' = H and R' 


R 

0 


(23) 


The material constants are 


%’ € o’ n ’ A m > R and H ( 24 ) 

The constant a Q is a reference stress that can be chosen conveniently in the 
stress range of interest. Note that for a « 0, i.e., near the virgin state of the 

material, eq. (21) gives e = e Q for a = a Q ; thus the constant e Q represents 
the initial creep rate at the reference stress a . 

ISOTROPIC LIMIT 

The isotropic limit of eqs. (11) — (15) is found by taking D-. to be an 

1 ^ 

isotropic tensor with = 1, i.e., D.. = ^ £... Under this limit, there results 


— _ 4 J 

7 s 2 ’ 

G - 75 4 

(25) 
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(26) 


Thus, the present creep model reduces to a Jg theory in the isotropic limit. 
Under uniaxial tension, the governing equations reduce to: 


'o 


(7— a 


^ 2n 


a o ' 


(27) 
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in which 



^=4: H a'^ 

st = 4- R 


We identify the behavior of the matrix material with that corresponding 
to the isotropic limit and observe that if we experimentally determine the 

constants o' , e^, n, 0, m, Si and <%! in eqs. (27) and (28), i.e., for the 

matrix, then using eqs. (29) we can find the constants o Q , c , n, (3, m, R and 
H in the multiaxial theory for the composite expressed in eqs. (11) - (15) and 
(20). Thus, we can characterize the multiaxial creep response of the composite 
knowing the uniaxial creep response of the matrix. As the effect of the fibers 
amounts essentially to a kinematic constraint, it is the matrix material that 
controls the inelastic behavior of the composite. 

In the following section, we outline a procedure for determining the 
constants for the matrix material. Following that, we make application of this 
procedure to Kanthal, a model matrix material of interest for high-temperature 
applications. 

DETERMINATION OF MATERIAL CONSTANTS FOR THE MATRIX 

Determination of the (isothermal) material constants o^, e' Q , n, 0, m, Si 
and <?€ from tests on the matrix material requires constant stress creep tests 
spanning the stress range of interest (fig. 2) and a sequence (at least one) of 
variable stress tests in the form of step-stress tests (fig. 3). Here, o' is a 
reference stress for the matrix material to be selected within the stress range 

of interest and is the initial creep rate (a « o) under o' Q (fig. 2). 

We note that the second of eqs. (29), with reference to the flow laws of 
eq. (21) and eq. (27), indicates that the ratio of the reference stress o' Q 
corresponding to the isotropic limit (identified here with the matrix material) 
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to the reference stress a of the composite under transverse tension is 



( 30 ) 


This result is anticipated in that the isotropic matrix material behaves as a 
( von-Mises ) material and with the addition of strong fibers is constrained to 
deform in transverse tension as a maximum shear stress ( Tresca ) material. 
The creep (or yield) strengths of these classical material idealizations are in 

the ratio fl/2. 

Now focusing on a typical step test (fig 3), we identify the creep rates 

immediately before and after the abrupt change in stress from to as 

and respectively. Assuming the stress change is made abruptly enough so 
that a remains essentially constant during the change, we obtain from eq. 27, 
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(31) 


Everything in eq. (31) is obtained from experiments except for the exponent n, 
therefore we can determine n. If data from several step tests are available, 
providing several independent expressions like eq. (31), we can find an optimal 
value of n satisfying these equations, e.g., in a least squares sense. At this 

stage, the flow law eq. (27) is fully determined as we know a' Q , e ' and n. 

Now we turn to the constant stress creep tests and focus on steady state 

creep (fig 2). The experimental data provides a steady state creep rate e for 

s 

each stress level a. The flow law, eq. (27), allows us to determine the steady 

state value a g of the internal state variable corresponding to each pair (c s , a), 
i.e., 


a 




e o 


1 

sir 


(32) 
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At steady state a = 0 and the evolutionary law eq. (28) gives 


e o 


St 2 m 
TF “s 


( 33 ) 


The only unknowns in eqs. (32) and (33) are m and the ratio St/ 76. As 
before, optimal values of these can be found satisfying several pairs of 
equations like (32) and (33) resulting from the available test data. We now 

know cr , e^, n, m and the ratio Stj 76. 

Now we focus on the primary creep stage. During early primary creep 
the first term in eq. (28) dominates and we can write 


(34) 

Upon integration of eq. (34) with the initial values £ = 0 and a = 0 we 
obtain: 


a = 


a 


76 



Again, from the flow law, eq. (27) 


(35) 


n (36) 

Now we consider a typical data point P during early primary creep as in fig. 

2. At P we know the creep strain e, the creep strain rate e , and the stress 
a. Thus for a given P the only unknowns in eqs. (35) and (36) are /? and 
76 . Again, consideration of a number of data points P during early primary 
creep yields optimal values of the constants /? and <76 . 



9 



This completes the procedure for specifying the material constants o' Qi e^, 
n, (3, m, 5£ and e# for the matrix material. As indicated earlier, use of eq. 
(29) then allows determination of the constants for the composite material, eq. 
(24). In the next section we present the results of applying this procedure to 
Kanthal, the matrix material of a Kanthal composite. 

DETERMINATION OF THE MATERIAL CONSTANTS FOR KANTHAL 

The Kanthal matrix material considered here is of the following 
composition by weight percentage: 73.2% Fe, 21% Cr, 5.8% A1 and 0.04% C. 

The data base for Kanthal is comprised of a set of constant stress creep tests 
(fig. 4) and a variable stress test (fig. 5) involving two abrupt stress steps. 
The data are isothermal at 650 C; this corresponds to T/T m s .45 for 
Kanthal. The constant stress tests are at 10,13,17 and 25 ksi (69,90,117 and 
173 MPa) taken over the relatively short time period < 2.2 hr with 
accumulated creep strain <1%, as seen in fig 4. The short time duration of 
the creep data exemplifies the cycle duration in some aerospace applications. 
The creep curves show typical primary creep and the minimum creep rate (the 
darkened portions of the curves) is identified as steady state for the present 
purposes. 

The step creep test of fig. 5 begins at 13 ksi (90 MPa) with a step to 
15.4 ksi (106 MPa) at 2 hr and a second step to 17.9 ksi (124 MPa) at 4 hr; 
the loading rate is 9000 ksi/hr (62100 MPa/hr). It is assumed that the stress 
steps are made sufficiently abruptly so that the negligible changes occur in the 
internal microstructure governing creep behavior (measured phenomenologically 
by a in eqs. (27) and (28) ). 

The test data was processed as outlined in the foregoing section. The 

reference stress o' Q was chosen as 10 ksi (69 MPa) with e ^ the measured 
initial creep rate at 10 ksi. Inelastic strain rates were estimated before and 
after each of the two stress steps of fig. 5 giving two equations, eq. (31), from 
which an optimal value of the exponent n was obtained. This was 

accomplished using a Levenberg/Marquardt least squares method. 

Steady state (minimum creep rate) data (f ,<r) were obtained from the 

s 

constant stress creep curves of fig. 4 and used in eqs. (32) and (33). Using 
the same least squares method indicated above the optimal values of m and 
55/ <%! were obtained. 
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Primary creep information was based on the three data points on the 10 
ksi (reference stress) creep curve indicated in fig. 4. Respective measurements 

of e, e and a at these points were used in eqs. (35) and (36) similarly 
providing optimal values of the constants /? and <%>. This completes the 
specification of the material constants for the Kanthal matrix. 

Values of the material constants for Kanthal and for a composite at 

650 °C are given in Table I (the material constants values are consistent with 
the units ksi and hr). 

Figs. 6 and 7 show correlations of the data of figs. 4 and 5 with 
calculations using the matrix constants of Table I in eqs. (27) - (29). Equiv- 
alently, the calculations could be made using the composite constants of Table 
I in the multiaxial forms eq. (11) - (15), (20) under the isotropic limit D-. -» 

K- 

PREDICTIONS OF THE CREEP THEORY 

The predictions of this section are made using the composite constants of 
Table I in eqs. (11) - (15), (20). Fig. 8 compares the creep response of the 
Kanthal matrix (isotropic limit) at 17 ksi (117 MPa) and 650 C with that 
predicted for transverse creep of the Kanthal composite, i.e., after the inclusion 
of strong unidirectional fibers (for example W-fibers) that suppress creep along 
their direction. As indicated earlier this amounts to a kinematic constraint 
that forces the matrix material, which would deform inelastically as a Jg 
( von-Mises ) material without fibers, to deform as a maximum shear stress - 
Tresca material. This is reflected in the ratio of reference stresses stated in 
eq. (30). Note, from eqs. (21), (27) and (29) the ratio of initial creep rates in 
fig. 8 is 


£ trans 
e matrix 

Here, with n = 3.3 this ratio is approximately 0.45. 

The lower curve (transverse creep) in fig. 8 can be interpreted as the 
creep response of the constrained matrix with no weakening from the 
introduction of fiber/matrix interfacial imperfections or from the evolution of a 
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weak interphase material through diffusion. In this sense, the lower curve of 
fig. 8 represents an upper bound on the transverse creep strength of the 
composite. Just as a comparison of creep of the unreinforced matrix with that 
of the composite along the fiber direction provides a measure of the 
effectiveness of fiber strengthening, a comparison of matrix creep with that 
transverse to the fibers (fig. 8) provides a measure of the integrity of the 
composite. Departure from the prediction of fig. 8 may indicate possible 
interfacial imperfection or degradation of properties resulting from the 
introduction of the fibers. 

A generic point in a structure composed of a unidirectional (or 
multidirectional) composite generally experiences stress components in addition 
to that aligned directly with the fibers; in fact, these adverse stress 
components (shear, transverse tension, etc.) control the inelastic response 
(creep) of the structure. This is illustrated by the final predictions shown in 
fig. 9. 

Fig. 9 shows the predicted responses of (closed end) thin walled Kanthal 
composite tubes subjected to an interior pressure. In one case the tube is 
axially reinforced (insert A in fig 9) and the other is circumferentially re- 
inforced (insert C). Each is subjected to an interior pressure p such that a = 
dR R 

jrf with y>>l. The axial stress is a z = a, the circumferential stress is Oq 
= 2 a and the radial stress is < 7 r » 0. 

For the axially reinforced tube (A), the governing eqs. (11) - (15), (20) 
become 


•A 


e 


z 


= 0 ; 




a = 


H 1 

5 ** 



R n m -0 
~ T G 



(38) 


Integration of these equations using the composite constants Table I gives the 
creep curve labelled €^(t) in fig. 9. 
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The corresponding equations for the circumferentially reinforced tube (C) 
are 



• r; 

These lead to the curve labelled e (t) in fig. 9. Clearly, the creep strength of 
the circumferentially reinforced Kanthal composite tube is much greater. The 
ratio of initial creep rates is 

= 2 21 x 97, (40) 

i>) 

differing by almost two orders of magnitude. This large difference reflects that 
the transverse stress controls the creep behavior. Transverse tension in the 
axially reinforced tube (A) is = 2<r ; that in the circumferentially 
reinforced tube (C) is a — a. In either case the fibers perform as expected, 
• A . c 

i.e., ensuring = 0 for case (A) and e q = 0 for case (C). 

SUMMARY AND CONCLUSIONS 

An anisotropic creep model is presented for metallic composites having 
strong fibers and low to moderate fiber volume percent (< 40%). 

Identification of the matrix behavior with that of the isotropic limit of the 
theory allows characterization of the composite through creep tests on the 
matrix alone. This amounts to assuming that the fibers provide a constraint 
that suppresses creep in their direction and, under transverse tension, forces 
the J 2 ( von-Mises ) matrix material to deform as a maximum shear stress 
( Tresca ) material. The fiber-matrix interface is assumed to be perfect without 
any weakening through the introduction of fiber/matrix interfacial imperfections 
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or from the diffusion related evolution of a weak interphase material. The 

prediction of transverse creep under these assumptions (fig. 8) thus provides an 
upper bound on the transverse creep strength of the composite. Comparison of 
the measured transverse creep response of the composite with that of the 
matrix as in fig. 8 provides an assessment of the integrity of the composite. 

A procedure is prescribed for characterizing the composite using creep 
tests on the matrix material. The required (isothermal) data base includes 

constant and step-wise constant creep tests. The procedure is applied to 

Kanthal composite, a model high-temperature composite material. Creep tests 

on Kanthal at 650 °C in the stress range 10-25 ksi (69-173 MPa) comprise 
the data base from which the material constants are derived. 

As an application of the theory, predictions are made of the creep 

response of thin-walled Kanthal composite tubes (closed ends) subjected to an 
interior pressure. Longitudinally and circumferentially reinforced tubes are 
considered. The predictions substantiate that circumferential reinforcement is 
far more efficient in limiting creep. The applications illustrate that creep is 
controlled, not by the stress along the fiber direction, but by the transverse 
component. This stress component is larger by a factor of two in the 
longitudinally reinforced tube. 

When weakening mechanisms arising from imperfect bonding, diffusion or 
the fiber volume percent becomes significant and the effect of the fibers is not 
limited just to a kinematic constraint as described earlier, the composite can- 
not be characterized accurately on matrix testing alone. The composite should 
then be modeled and tested as a material in its own right as in [4,5]. Then 
the intrinsic weakening effects and their time and history dependent evolution 
are appropriately reflected in the experimental results. 
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Table I. Material Constants 


Matrix 
= 10 

e' = .00165 
n = 3.3 
0 = 1.05 
m = 5.78 

9t = 9.5 x 10 
<#= 38.8 


Composite 

cr = 11.547 
o 

e n = 0.00143 
o 

n = 3.3 
0 = 1.05 
m = 5.78 

R = 3163.68 
H = 0.3558 
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(b) TRANSVERSE STRESS. 

FIGURE 1. - MATERIAL ELEMENTS UNDER STRESS. 



TIME 



FIGURE 2. - TYPICAL CONSTANT STRESS CREEP CURVES FOR THE MATRIX 
MATERIAL. Og' IS A REFERENCE STRESS AND££ IS THE INITIAL CREEP 
RATE AT O J 0 . P IS A GENERIC POINT IN EARLY PRIMARY CREEP. 



a 2 

°1 



TIME 



FIGURE 3. - TYPICAL STEP-STRESS CURVE FOR THE MATRIX 
MATERIAL. 6, IS THE FINAL CREEP RATE AT 0, (BEFORE 
STEP), C 2 IS THE INITIAL CREEP RATE AT 0 2 (AFTER 
STEP). 



TIME, HR 

FIGURE N. - EXPERIMENTAL CREEP CURVES FOR KANTHAL AT 10, 13, 
17. AND 25 KSI AT TEMPERATURE 650 °C. MINIMUM CREEP RATES 
(DARKENED PORTIONS) ARE TAKEN AS STEADY STATE. PRIMARY 
CREEP INFORMATION IS BASED ON THE THREE POINTS INDICATED 
ON THE 10 KSI (REFERENCE STRESS) CREEP CURVE. 
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TOTAL STRAIN, PERCENT CREEP STRAIN, PERCENT TOTAL STRAIN, PERCENT 



0 12 3 4 5 6 


TIME, HR 

FIGURE 5, - EXPERIMENTAL STEP-STRESS CREEP CURVES FOR KANTHAL 
AT 13-15.4 KS I (40-106 MPa) AND 15.4-17.9 ksj (106-124 MPa) 
CREEP RATES ARE IDENTIFIED BEFORE AND AFTER EACH STEP. 



TIME, HR 


FIGURE 6. - CORRELATIONS OF EXPERIMENTAL (CIRCLES) AND CAL- 
CULATED (LINES) CONSTANT STRESS CREEP CURVES. FOR 10, 13, 
17, AND 25 ksi (69, 90, 117, AND 173 MPa) AT 650 °C. 


10 KSI (64 MPa) 



0 123456789 


TIME, HR 

FIGURE 7. - CORRELATIONS OF EXPERIMENTAL (CIRCLES) AND CAL- 
CULATED (LINES) STEP-WISE CREEP TESTS AT 13-15.4 ksi (40- 
106 MPa) AND 15.4-17.9 ksi (106-124 MPa). ALSO SHOWN IS 
A PREDICTION OF A FINAL STEP 17.9-10 KSI (124-64 MPa) AT 
650 °C. 



0 0.5 1.0 1.5 2.0 


TIME, hr 

FIGURE 8. - COMPAR I SON OF PREDICTED MATRIX CREEP AND TRANS- 
VERSE CREEP OF THE COMPOSITE (UPPER BOUND ON TRANSVERSE 
CREEP STRENGTH) AT 17 KSI (117 MPa) AND 650 °C. 



0 0.5 1.0 1.5 2.0 

TIME, hr 

FIGURE 9. - COMPARISON OF PREDICTED CREEP RESPONSES FOR 


AXIALLY REINFORCED (A) AND CIRCUMFERENTIALLY REINFORCED 
(C) KANTHAL COMPOSITE TUBES UNDER INTERIOR PRESSURE 
CORRESPONDING TO 0 = PR/2t = 7 KSI (48.3 MPa). 
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